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Abstract 

A phenomenological approach to the ferromagnetic two dimensional 
Potts model on square lattice is proposed. Our goal is to present a simple 
functional form that obeys the known properties possessed by the free 
energy of the q-state Potts model. The duality symmetry of the 2D Potts 
model together with the known results on its critical exponent a allow to 
fix consistently the details of the proposed expression for the free energy. 
The agreement of the analytic ansatz with numerical data in the q — 3 
case is very good at high and low temperatures as well as at the critical 
point. It is shown that the q > 4 cases naturally fit into the same scheme 
and that one should also expect a good agreement with numerical data. 
The limiting q = 4 case is shortly discussed. 
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Introduction 

The Potts modeQ is the most natural generalisation of the two dimensional Ising 
model and it is deeply connected with many areas of both statistical physics and 
mathematics (for nice detailed reviews see [I] , [2] , [3] , [3] and [5] ; for the relation 
with the problem of colour confinement see [Hj and the review [7]). 

It has been very difficult at the moment to compute directly the analytic 
free energy of the ferromagnetic Potts model in two dimensions on square lattice 
(which is the main object of the present paper) as it was done by Onsager in 
the Ising case [5]. Thus, it will be proposed here a less direct approach which 
is, nevertheless, able to give non-trivial analytic information on the free energy. 
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x In the present paper we will only consider spin models on regular square lattices. 
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The main purpose of this approach is not to give an exact solution for the free 
energy, but rather to present a simple functional form that obeys the known 
properties possessed by the free energy of the q-state Potts model. In particular 
this method may give a reasonable approximation of the free energy function 
in all the range of temperatures even the not yet explored ones, i.e. outside 
the results of high and low temperature or the neighbourhood of critical point. 
Following the point of view of the phenomenological Regge theory of scattering 
(see [3] and [TU]; two detailed reviews are [TT] and [H]), a reasonable analytic 
form for the free energy in terms of one free parameter will be derived. In the 
present case it plays the role of Regge's trajectories in high energy physics since 
it parametrises some analytical features of the Potts model such as its duality 
properties and the Fisher zeros, in an analogue way as the Regge trajectories 
encode the non-perturbative duality properties of scattering amplitudes, as was 
first observed in [13] . This can be done by requiring that the sought analytic 
expression for the free energy of the 2D Potts model should be compatible (in 
a suitable sense explained in the next sections) with the proposal made in |14j , 
[T5] , [TB] for the free energy of the three dimensional Ising model. The proposal 
of [14], [15], [16] is in a very good agreement with numerical data so that one 
should expect that this approach should provide one with a formula for the free 
energy of 2D Potts which fit equally well the available numerical data. It will 
be shown that this is indeed the case. A non-trivial by-product of the present 
approach is that the proposed free energy of the 2D Potts model has the locus of 
the Fisher zeros which coincides with the well known and well tested conjecture 
in [17]. 

The paper is organised as follows: in the second section, it will be discussed 
a suitable consistency condition (called here " dimensional compatibility" ) which 
allows to derive an analytic ansatz for the free energy of the 2D Potts model in 
terms of few q— dependent curves. In the third section, it will be analysed how 
the known duality symmetry of the Potts model and its known critical behaviour 
fix all but one curve. In the fourth section, the proposed ansatz will be compared 
with the available numerical data both at high and low temperatures as well 
as at the critical point in the case q — 3. In the fifth section, the q > 4 cases 
will be discussed. In the sixth section it will be shown how the proposed ansatz 
automatically predicts the locus of Fisher zeros consistent with the well known 
conjecture. 

1 Dimensional Compatibility 

It is well known that the main thermodynamical quantities are combinatorial 
in nature. This fact makes the analysis of three-dimensional lattices even more 
difficult than the two dimensional cases. In this respect, it is a rather surprising 
result that one can even imagine to write down a simple explicit functional form 
for the free energy of the Ising model in three dimensions which obeys known 
physical properties (in particular, good agreement with the available numerical 
data). Such simple expression is based on the assumption that the change in the 
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combinatorial complexity when going from one to two dimensions is formally 
quite similar to the analogous change in the combinatorial complexity when 
going from two to three dimensions. One can write the exact Onsager solution 
for the free energy in two dimensions as a convolution integral of the known 
exact free energy in one dimension with a suitable kernel. This integral kernel, 
which can be found explicitly, is the mathematical object responsible for the 
change in the combinatorial complexity when going from one to two dimensions 
[14] . Then, one argues that the change of many combinatorial quantities on 
hypercubic lattices (such as, for instance, the scaling of the number of states of a 
given energy with the size of the lattice) when going from one to two dimensions 
are similar to the analogous change when going from two to three dimensions 
(this is somehow confirmed by a comparison with the available numerical data 
[16] , [15] )■ Thus, it is not unnatural to assume that a suitable integral kernel 
exists such that when one consider its convolution with the Onsager solution one 
gets the (still unknown) exact free energy of the Ising model in three dimensions. 
Our assumption based on the above considerations is that such kernel is as 
similar as possible to the already known one allowing to jump from one to two 
dimensions. 

The Kallen-Lehmann representation [T¥ gives rise to an ansatz for the free 
energy of the 3DI model of the following for 
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where the values of the parameters in the case of the two-dimensional Ising 
model would be A = 1, Ci = 1/2 for i = 0, .., 3. 

Namely, one may find the operator Di^ (see Eq. ([2])) which "dresses" the 
trivial one-dimensional solution of the Ising model giving rise to the Onsager 
solution. Then, the Kallen-Lehmann free-energy for the three dimensional Ising 
model is deduced [14] by modifying D 12 in such a way that the parameters 
which in Di 2 are fixed to be 1 and 1/2 (namely, in the two dimensional case, 
A = 1, ( { = 1/2) become the free parameters. Afterwards, one can use such a 
modified operator to dress the Onsager solution obtaining a useful ansatz for 
the free energy of the three dimensional Ising model. To be more precise, it is 
possible to define a class of operators T>D.q such that when they act on the free 
energy Fd q of the q— state ferromagnetic Potts model on hypercubic lattices in 



2 The parameters appearing in these formulae are related with those appearing in |14| by 
the following identities = u, f 3 = 1/2 = f 2 , Co = a - 
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D dimensions give rise to the free energy Fm+i),q of the q— state Potts model 
in (D + 1) dimensions 

Dd,, 08)] = F (JJ+1) , 9 08) • (2) 

To fix the arbitrariness of the above definition it necessary to specify the domain 
of T>D.q and how it acts. Since we want that the free energies of reasonable 
systems belong to the domain of T>D, q , the domain of Dd,, will be the class 
of functions which are smooth on R + besides at most a finite number of points 
in which some of the derivatives of the function may be discontinuous. The 
variable is the inverse temperature and the singular point (or points) represents 
the phase transition. The range of the operator coincides with the domain since, 
by definition, when this operator acts on the free energy of the Ising model on 
a hypercubic lattice it generates the free energy of an analogous system in one 
more dimension. A reasonable way to represent these operators is as integral 
Kernel [H] : the simplest of these integral kernel can be constructed explicitly 
by comparing the Onsager solution with the free energy of the one-dimensional 
Ising model. Then, one can use an integral Kernel of the same functional form to 
act on the Onsager solution obtaining an approximate ansatz for the functional 
form of the free energy of the Ising model in three dimensions (the comparison 
of the ansatz with the numerical data is quite promising [15j |16j). 



1.1 Consistency condition 

Let Qd,9 be the operators which when applied to the free energy Fu, q of the 
q— state ferromagnetic Potts model in D dimensions give rise to the free energy 
Fr>,( q +i) of the (q + 1) —state Potts model in D dimensions 

Qo, q [F D , q 08)] = F D>(q+1) 08) . (3) 

The domain of Q^.g is the class of functions which are smooth on R + besides at 
most a finite number of points in which some of the derivatives of the function 
may be discontinuous (again the singular points represent the phase transitions 
of the system). Also in this case the range of the operator coincides with the 
domain. Therefore, being the ranges and the domains of the operators T>d q 
and Qd. 9 compatible, it makes sense to compose them. In particular, one can 
observe that they have to satisfy a sort of commutativity constraint: 

Q(Z»+l),g " T*D, q = D Di ( g+ i) • Qz),g- (4) 

If one supposes that Eq. ([T]) is the exact free energy of the three dimensional 
Ising model for some values of the parameters, the simplest way to satisfy Eq. 
Q is the following approximate analytic formula for the free energy of the two 
dimensional ferromagnetic Potts model: 

F 2D (q,u) = C q + ^-J dt log 



(k 2D (q,u)) 2 < 1; C g = log ^ 2 ex P(^) + (g- 1 ) ^ . u = exp ( _ /8)) (5) 



/ 9 „ \ ?(<?) 
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where km (<?, u) is the function encoding the duality properties of the model, 
the function C q can be found by comparing the high temperatures expansion of 
Ising and Potts models in two dimensions (see, for instance, [I]), A (q), <, (q) are 
two parameters 0. It is worth to note that we are still free to add a constant 
(which depends on q) to the above free energy. Unfortunately this fact prevents 
us from using the Baxter's results on the free energy on the critical point to fix 
X(q). It is trivial to verify that the critical point u c corresponding to the free 
energy in Eq. ([5]) is determined by the equation 



as one expects on the basis of the Onsager solution. 

When q < 4, the curve ^ (q) and &2.d(<3S u) can be fixed a priori using theo- 
retical arguments related to the duality symmetry and to the known results on 
the critical behaviour. While we will fix the normalisation X(q) by a comparison 
with the numerical data. 

As it will be shown in the next sections, ^(q) is related to the critical exponent 
a (q) . Through the well known critical behaviour of the Potts model, one can 
get an implicit functional relation between c; and a 



In the q > 4 cases (in which the transition is first order) one can fix a priori 
<; (q) in terms of A (q) using the known exact results of Baxter on the latent heat. 
Indeed, the above formula may look an ad hoc approximation, or, at least, not 
very natural at a first glance. On the other hand from the arguments at the 
beginning of section[T]it stems a constraint on the free energy of the Potts model 
in two dimensions. Namely, we are interested in finding a functional form that 
exactly obeys certain properties known to be possessed by the true free energy. 
So, we search for an expression for the free energy which is compatible with 
the properties of the integral kernels discussed above. This gives a further 
constraint on the form of the free energy of the Potts model in two dimensions. 
From the combinatorial point of view, this makes the proposed form in Eq. ([5]) 
natural. In other words, besides the known constraints on the free energy of 
the two dimensional Potts model (such as the critical exponent and the duality 
symmetry), the proposed form in Eq. ([5]) is also the simplest compatible with 
the recursive structure proposed in [14] (which has been proved to be in good 
agreement with known results [16] , |15j). This explains why it is quite useful 
to look at the recursive structure connecting the Ising models in two and in 
three-dimensional to obtain an additional bit of information in order to fix the 
residual arbitrariness left in the choice of the free energy. 

3 X(q) is related to the overall normalisation of the non-analytic part of the partition function 
while <; (g) is related to the critical behaviour. 



k2D(q,U c ) = 1 



(6) 



C (?) = $(*(«)). 
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2 Duality and /^(^w) 



The duality transformation in the case of the two-dimensional Ising model, was 
discovered in [IS] before the exact solution of Onsager [8], In the case of the 
Potts model we are considering the following Hamiltonian: 

H = -J <W; > ( 8 ) 

<ij> 

and the duality is 

D(u)= )~ U . (9) 
1 + (q - 1) u 

Note that this transformation is not a symmetry of the full free energy per site 
in the thermodynamic limit: it leaves invariant the non- analytic part of the free 
energy while the trivial term log(- — ) is not invariant. Anyway the critical 
point is determined by the properties of the non-analytic part. The fixed point 
of the duality transformation u c — u c (g) is 

D (tt e ) = u c u c (q) = 1 (10) 
i + v<7 

Thus, one has to find a function k2D(q,u) which encodes the duality properties 
of the two dimensional Potts model and which reduces to the known result when 
q = 2: that is 

k 2 D(q,D(u)) = ±k 2D (q,u) (11) 

where the ± signs appear because the Kallen-Lehmann free energy in Eq. §5§ 
depends on (&2.d(<?, P)) ■ 

The simplest solution (let us call it k) of Eq. ([TT|) is 



it 



{\-{q-l)u 2 + {q~2)u) 



k(q,u)=A —— 2 '-. (12) 

(1 + {q — l)u z ) 

it can be easily seen that k(q, u) fulfills Eq. (fTT|) for any value of the constant 
A which we will fix with the normalisation condition at the critical point in Eq. 

m 

k(q,u c ) = 1 =*> 

1 _ U C (1 - (q - 1) U 2 C + (g - 2) M C ) 

^(«) " (l + ( 9 -l) M 2) 2 

Furthermore the critical point is located at u — u c in Eq. (fTCTj) as it should be 
and it can be easily seen that when q = 2 it reduces to the expression of k 2 D in 
terms of the low temperatures variable u. 

Indeed, once one has found the simplest k invariant under duality transfor- 
mation and which reduces to the Ising case when q = 2, one can construct many 
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more solutions by simply taking functions / of the k in Eq. (fT2"j) 

JM* «) = / U (?) M(1 ' (g ' 1)u2 + (g 2 ' 2)u) l 

V (l + (g-l) U 2 ) 2 ; 

such that / (x) has the maximum when x = 1 and /(l) = 1. The simplest 
possibility is to consider / (x) of the form 

/ [x) = x E ^ 

One procedure to determine E (q) for q — 3 is to look at the coefficients of 
the low temperatures expansions (see [20], [H])- Using Eq. one can compute 
the ratios S n for very small u 

c / \ a n (?) ") 

o n {q,u) = 



a n (q, u) = 



a n+1 (q, u) ' 
d n F 2D (q,u) 



du 7 ' 



for some small n and verify that it is possible to fulfil the expected scaling in 
the cases q = 3 with the choice 

E(3) = 2, (13) 

so that we will take 

u (l - 2u 2 + u)\ 



k 2D (3,u)= A (3) a , 

y (1 + 2m^j y 

while, of course, the compatibility with the Onsager solution tells that E (2) = 1. 
It is also interesting to observe that consistency with the q — » 1 + limit (where 
the u-dependence disappears) would suggest E(l) = 0. 
Thus, from now on, we will fix E (3) as in Eq. ([13")) . 

2.1 The critical behaviour 

The critical exponent a (q) for two dimensional Potts model (when q < 4) is 
known to be (see [T5] ) 

2(1 -2a) 
3 (1 — x) 

x = — arccos - — 
7T V 2 

where the positive values of x correspond to the tricritical point while the neg- 
ative values correspond to the critical point. 
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One can fix a priori the curve ^ (q) by looking at the critical behaviour of 
the model: the specific heat is known to have (see [22], [23]) the following forms 
for q = 3 

C dm {q = 3,u^ u c ) <x\u-u c (3)P 1/3 . (14) 

On the other hand, the second derivative of the (non-analytic part of the) free 
energy F2D in Eq. ([5]) reads 

02 p _ *(g)A(g) /oa^ / (sin 2 x)A^)-^ 



2vr K ' J 1 + A«W V 1 + 

where 



A = l-ffsin 2 ^ 0<c(g)<l, 
# = (feD(u)) 2 . 

Near the critical point, if « 1 and (9 u if « since, as it has been already 
discussed, the critical point is a smooth maximum of k,2D (u). For this reason, 
the most singular term is the first one: 



2ir v u > J 1 + A ? («) 


since the divergent term in the integral of the second term in Eq. (II 5|) are 
compensated by the vanishing first derivative of H at the critical point (while 
d^H is of order 1). By imposing that the singular part of the specific heat in Eq. 
(|16p reproduces the known critical behaviour in Eq. (|14[) one gets an implicit 
relation between a and ? which fixes <; (3) to be 

?(3) = 0.4 (17) 

In the next sections we will draw a picture of the critical part of the free energy 
in Eq. ([5]) against the known results at the critical point which shows a excellent 
agreement. 



3 Comparison with numerical data for q = 3 

An explicit analytic expression for the free energy in Eq. ([5]) of the two di- 
mensional ferromagnetic Potts model has been constructed in which theoretical 
arguments (basically, duality and the known critical behaviour) can fix every- 
thing but one curve X(q). Indeed, it is easy to see that A (3) can fixed in 
terms of the numerical expansion data at low temperatures giving rise to a 
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good agreement. We would like to further emphasise that the main purpose of 
our work is to obtain an unique approximate functional form for free energy in 
the widest possible range of temperatures. This is conceptually different and 
complementary point of view with respect to the known numerical analysis and 
perturbative expansions. So we have not just assumed the soundness of the low 
temperature and the critical point numerical studies, but we used them to fix 
the only free parameter of the ansatz. As it is explained below the precision 
of our numerical algorithm is not sensitive as the numerical data themselves, 
neither we pretend to compete with the precision of these methods. Rather 
we mean to give a global description for the behaviour of the free energy in a 
unified scheme compatible with the known results on their respective different 
domains of applicability. Despite the functional simplicity of the free energy in 
Esq. ([5]) a low temperature Taylor expansion, which involves many numerical 
derivations and integrations, using the standard commercial software (available 
to us) is very inaccurate and anyway beyond the goal of this work. 

3.1 The low and high temperatures 

Because of the built-in duality invariance of the model, we will only need to check 
the agreement at low temperatures since the agreement at high follows from 
duality. Our analysis is based on the references [3D], [3T]: to be more precise, 
we checked that the normalisation in |21) is consistent with (and reduces to) 
the normalisation of [3U] m the q = 2 case of the two dimensional Ising model. 
In particular, the expansion of the partition function in the above references 
corresponds to only consider the interesting non-analytic term (neglecting the 
trivial term log cosh 2(3 in the Ising case). Therefore, in the Potts case one has 
to compare the low temperature expansions for q = 3 with the second term on 
the right hand side of Eq. (|5|). We consider the low temperature expansion, 
up the 14 th order, valid no more than u c (3)/100, since already in the 2D Ising 
model the Onsager solution deviates significantly from the numerical expansion 
when u > u c (2)/100. The low temperatures expansion of the free energy in [20] . 
[2Tj for the q = 3 case is 



For this value of A (3) at low (and therefore high) temperatures the "precision" of 
the Kallen-Lehmann ansatz versus the numerical free energy may be measured 
in many different ways; for instance one can use the following: 



F mc {u) 



Log(l + 2u A + 4u 6 + 4u 7 + 6u 8 + 2Au 9 + 

+24u 10 + 68U 11 + 190u 12 + 192u 13 + 904u 14 ) 



and the best value we have found for A (3) is 

A(3) = 0.1543 




As one can see from Fig. (|3.1[) the agreement is very good. 
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Figure 1: Numerical F nm (u) and Kallen-Lehmann i«2£j(u) free energies at low tem- 
peratures for g = 3. 

3.2 The critical behaviour 

Because of the simplicity of the free energy in Eq. ([5]). one can perform the 
numerical graph using Mathematica® . The picture |2] , in which we plot the 
Kallen-Lehmann free energy F2d(3,u) against the critical free energy 

Fcritic = a + b \u - u c (3) | 3 

in a neighbourhood (sized 1% of u c (3)) of the critical point show a remark- 
able agreement. Also in this region the precision p(F cr i t i C , F^d) of the result 
maintains well below 1%. 

To get a clearer idea on the precision of the proposed ansatz, one could 
compare, using Mathematica®, the exact Onsager solution at high temperature 
versus the corresponding high temperatures expansion in |20j in an interval of 
u from to u c (2)/100 (the Onsager solution deviates significantly from the 
numerical small u expansion when u > « c (2)/100). If one would do this, one 
would recognise how well the proposed ansatz describe the available numerical 
data in the 2D Potts case for q = 3. 

4 The q > 4 cases 

The q > 4 cases are qualitative very different from the cases q < 4 since, when 
q > 4, the phase transition is of first order. Thus, one could expect that even 
if the ansatz in Eq. ([S]) works very well for q < 4, it is not at all obvious if it 
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Figure 2: Critical behaviour of free energy F cr uic and Kallen-Lehmann F2d(u) near 
the critical point (q = 3, a = -0.0147232, b = 0.677394) 



can also work when q > 4. Remarkably enough, it can be shown that the free 
energy in Eq. ([5]) does indeed describe first order phase transition at the critical 
point u c (q) in Eq. (|TU1) provided ? (q) becomes negative: 

?(<z)<o. 

When c (q) is negative (?(<?) = — |? (<z)|)j the leading behaviour of the internal 
energy corresponding to the free energy in Eq. ([5]) is E^d (q, u ~ m c ) « 



A(£) 
2tt 



(2k 2D d u k 2D ) 



(q) (l - (k 2D (q,u)) 2 sin 2 t 



-k(8)l-i 



sin 2 i 



(ft 



l+(l-(fc 2D (q,u)) 2 sin 2 i 



2tt 



7T 

(2k 2D d u k 2D ) J q(q) (l- (k 2D (q,u)) 2 sin 2 * sin 2 t 



dt. (18) 



Indeed, the derivative of k 2 u(q,u) at the critical point is still zero but the 
divergence of the integral is stronger so that it gives rise to a finite contribution, 
for u~ u c one has: 

\d u k 2D \ « 2e -C 1, fc 2 _D « 1 - e 2 , 
while the integral in Eq. (fT5|) is dominated by the region t tt/2. Thus, one 
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gets a finite discontinuity of the first derivative 

— a 

oo 

^ i i iy = T2X(qy{q) 

7r J i + y z 

— oo 

where a is a small positive number, the — sign refers to taking the limit to the 
critical point from the left and the + from the right. Therefore the internal 
energy acquires a finite jump (namely, a non-trivial latent heat) at the phase 
transition as it should be. Furthermore, the jump is symmetric in agreement 
with the result of Baxter. With the same arguments one can see that when q (q) 
is positive and less than one (as it is the case for q < 4), both left and right 
derivatives are zero at the critical point so that the transition is second order. 
This is a very interesting fact in itself since it allows to describe with the same 
analytical ansatz also the region in which a first order transition takes place. 

The above reasoning tells that to describe the first order region s (q) < 
and to describe the second order region <; (q) > 0. In the second order region it 
has been shown that ^ (q) is related to the critical exponent a. Below it will be 
discussed that in the first order region one can fix the discontinuity of the first 
derivative of the free energy at the critical point (which, roughly, is proportional 
to \(q)s (q)) using some exact results obtained by Baxter. 



1 - (1 - e 2 ) (1 - <T 



dS 



4.1 Comparison with the Baxter results 

Baxter (see [2], for a review see Q]) was able to compute the latent heat L (q, u c ) 
at the critical point in the q > 4 case: 

L(q,u c ) = 2 f 1 + q- 1 / 2 ) A (q) tanh ( | ), q > 4, 



A(q) = Yl (tanh(n6)) 2 ; cosh 9 = ^ 

n=l 



Thus, one can fix the discontinuity AE C of the first derivative of the free energy 
in Eq. (O at the critical point 

AE C = E 2 d (<?, w+) - E 2 d (q, u~) , 
in terms of the Baxter result: 



AE C = L (q, u c ) . (19) 

Such equation allows, in principle, to fix one of the two curves in terms of the 
other (for instance, one can choose to express A (q) in terms of <j (q)) living us 
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with only one curve which could be fixed by looking at the numerical expansion 
at low temperatures. Unfortunately at least in the cases q = 5 and q — 6 (which 
we analysed more closely) , we have not been able to develop a suitable software 
to test at the same time the low temperatures and the critical behaviour. Even 
if, at first glance, the numerical problems to be solved for q > 4 are similar 
to the ones which appear for q < 4, there are two important differences. The 
first is that for q < 4 the known critical behaviour is only related to <; (q) which 
therefore can be fixed, while for q > 4 the Baxter result on the latent heat 
determines AE C which is a rather complicated function of both A (q) and ? (q) 
and this does make the numerical analysis more involved. The second is that, 
at least in the cases q = 5 and q = 6, the numbers which arise in the low 
temperatures expansions are very small and this makes our software extremely 
slow. 

Nevertheless, we have verified using the software Mathematica® that look- 
ing at the low temperatures only, it is possible to achieves an almost perfect 
agreement with the series expansion of [5D] and [3T|, both for q = 5 and q = 6. 
This is a strong indication that this framework also works in the case q > 4 
since, because of the built in invariance under the duality transformation in 
Eq. ([5]), an excellent matching at low temperatures by construction implies an 
equivalent agreement at high temperatures. Thus, taking into account that in 
the s(q) < region the free energy in Eq. ([5]) has a first order phase transi- 
tion, one should expect that the present method provides an explicit analytic 
description of the free energy both for q < 4 and for q > 4 in terms of only two 
q— dependent curvefl in excellent agreement with numerical expansion data. 



4.2 The q = 4 case 

The q = 4 case is the more delicate. The first obvious reason is that q = 4 
is the boundary between the range in which the model exhibit a second order 
phase transition (q < 4) and the range in which the transition is first order 
(that is, q > 4). As a matter of fact, for q = 4: the specific heat singularity has 
logarithmic correction (see [22] or [23]): 

C div (q = 4,u- u c ) « |u " Mc(4)l ^ . (20) 
log |w — u c (4)| 

The present formalism provides one with a very natural mechanism for the 
arising of logarithmic corrections. 

Assuming that c; (q) is a continuous function of q, then one could argue that 
S (4) = since in the first order region q (q) < while in the second order region 
< ? (q) < 1. The interpretation of the last sentence is that, when q (q) = 0, the 
present framework provides one with the appearance of the expected "nested" 



furthermore, at least one of these two curves can be fixed a priori analytically using 
known analytical results at the critical point (see Eqs. 117t and ([19}). 
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logarithm inside the free energy: 



F 2D (4,u) = C i + 



A(4) 
2vr 



dt log 



l r 

2 



1 + S (q, u) log (l - (fc 2D (4, u)f sin 2 i) 



Unfortunately, we have still not found a theoretical argument to fix a priori E in 
an analytic way. Nevertheless, it is worth to stress that the above formula does 
give rise automatically to a second order phase transition with a logarithmic 
correction of the type in Eq. (|2"0)) and, therefore, to find theoretical arguments 
able to fix S (q, u) is an interesting open problem. 



5 Fisher zeros 

A powerful theoretical tool is the analysis of the Fisher zeros [23] in which the 
inverse temperature /3 is extended to the whole complex /3-plane (in the same 
way as Yang and Lee complexified the magnetic field [25 ). By looking at the 
distribution of zeros of the partition function in the complex /3 plane, one can 
determine the universal amplitude ratio A + / 04_ of the specific heat and write 
down simple expressions (which only involve the density of zeros and the angle 
which the line of zeros form with the real /3 axis at the critical point) for the 
free energy and the specific heat close to the critical point; see [UJ [55]. Such 
tools are also useful when analysing the strength of the phase transitions (see 
for instance [28], [29], [30] and references therein). Therefore, the Fisher zeros 
contain very deep non-perturbative informations on the corresponding systems. 

It is expected that also in the case of the two dimensional Potts model 
the Fisher zeros should lie on a circle: strong theoretical as well as numerical 
evidences have been provided in [17] (for some more recent evidences see [31] . 
[32] and references therein). The circle is given by \x\ = 1, where x = v/^Vq and 
v = eP — 1 = vT 1 — 1. To study the locus of the Fisher zeros, it is convenient to 
express k(q,u) as follows: 

~, . uD (u) . . 

%,n) = 4 \; a ■ (21) 
(u + D (uj) 

It is easy to observe from ([21]) thalU k(q, u) < 1 and k(q, u) = 1 iff u = D(u); this 
condition identifies the real fixed point of the duality map. On the other hand 
according to the conjecture, the locus of Fisher zeros of the 2D Potts model on 
square lattice is the analytic extension of the equation u c = D(u c ). Therefore 
our proposal is consistent with the conjecture providing further support to this 
framework. 



5 One may notice that in the above expression in Eq. I|21|l for k(q,u) the dependence 
on q is implicit in the duality transformation D(u). Furthermore, k(q,u) reduces to the 
known expression for q = 2. Therefore, when expressed in terms of u and D(u), the equation 
determining the locus of Fisher zeros is formally exactly the same as in the Ising case provided 
one replaces the duality transformation of the Ising case with the corresponding Potts duality 
transformation. 
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6 Conclusions and perspectives 



A phcnomenological approach to the ferromagnetic two dimensional Potts model 
on square lattice has been developed. After introducing the D— dressing and the 
q— dressing operators DD,q and Qd. 9 , it has been described how the compatibil- 
ity between T)r> _ q and Qd. 9 allows one to write down an explicit analytic ansatz 
for the free energy in terms of one free parameter (for each q). The duality 
symmetries of the 2D Potts model together with the known theoretical results 
on its critical exponent allow to fix a priori all but one curve. The agreement 
of the proposed analytic free energy with low and high temperatures expansion 
as well as the critical point is excellent for q = 3. For q = 5 and q = 6 one can 
also see that the agreement with numerical data at low and high temperatures 
is also very good but we have not been able to test the corresponding critical 
points because of some subtle numerical problems. Nevertheless it has been 
proved that the corresponding phase transition when <; < is hrst order. The 
q = 4 case remains basically opens but we have some indications that the present 
framework is also able to capture important features of such subtle case since 
it predicts automatically logarithmic correction to the power law divergence of 
the specific heat. 
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